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ABSTRACT 

In order to investigate how the growth of galactic bulges progresses with the 
growth of central black holes (BHs), we observed molecular gas (fuel for the 
coming star formation) in possibly young active galaxies, narrow-line Seyfert 1 
galaxies (NLSls). We present the results of radio observations of 12 CO (1 — > 0) 
using the Nobeyama Millimeter Array (with 2-4 kpc spatial resolution) for two 
FIR-bright NLSls, yielding the first detection of their CO emission. Correspond- 
ing molecular-gas masses M^Hq) of (1 — 3) x 1O 9 M are the 2nd and 4th largest 
ones among NLSls. By estimating dynamical masses and bulge masses Mb u i ge for 
these two NLSls using CO channel map and CO line widths, we found M(H 2 ) 
amount to 0.13-0.35 of these masses. Taking account the star formation efficiency 
(~ 0.1), the increase in Mb u i ge in those NLSls in the near future ( ^ 10 75 yr) is 
expected not to be a huge fraction (1-5 % of the preexisting stars). Bulge growth 
may have finished before BH growth, or bulge-BH coevolution may proceed with 
many, occasional discrete events, where one coevolution event produces only a 
small amount of mass growth of BHs and of bulges. We also discuss the ratios of 
star-formation rate-to-gas accretion rate onto BHs, finding that two NLSls have 
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very small ratios (~ 1) compared with the Mb u i ge /MBH ratios found in active and 
inactive galaxies (~ 700). This huge difference suggests either the non-overlapped 
coevolution, long star formation duration or temporarily high accretion rate dur- 
ing NLS1 phase. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: individual 
(IRAS 04312+4008, IRAS 05262+4432) - ISM: molecules - radio lines: ISM 

1. INTRODUCTION 

Mass of galactic spheroid component (galactic bulge or elliptical galaxy itself) is tightly 
correlated with mass of central black hole (BH) both in normal and active galaxies (Gebhardt 
et al. 2000; Ferrarese & Merritt 2000; Nelson 2000). Along with the lack of its cosmic 
evolution (Shields et al. 2003; Kiuchi et al. 2006; see however, Akiyama 2005, Woo et al. 
2006, and Borys et al. 2005), the tightness of the correlation indicates that the growth of 
spheroid mass and that of BH mass (M B h) occur almost contemporarily. 

However, no direct evidence for simultaneous coevolution has been found, and thus its 
physical origin remains a great mystery. Jet and/or outflow from an Active Galactic Nucleus 
(AGN) may govern the formation of spheroid (Silk & Rees 1998); i.e., BHs made before the 
bulges' growth. On the contrary, radiation drag by massive stars during the bulge formation 
may have controlled the BH growth (Umemura 2001), resulting in bulges formed prior to 
BHs' growth. 

Most numerical studies presume that cool, interstellar gas drives coevolution (e.g., Kauff- 
mann & Haehnelt 2000). For instance, semi-analytical approaches assume a major merging 
of galaxies initiates the growth of BHs and bulges, and usually adopt an AGN duration 
(e-folding accretion time scale) around 10 7 ' 5 yr and a star formation duration of ~ 10 8 yr 
(Kauffmann, White & Guiderdoni 1993; Kauffmann & Haehnelt 2000). Theoretically and 
observationally, it is unknown which growth phase precede another. Now, let us simply as- 
sume that the BH and bulge growth phases begin conjunctionally and proceed with similar 
durations. Then, we will see enough molecular gas - fuel for the coming star formation - in 
the bulges of AGNs in the early stage of BH growth (e.g., in the first quarter of the AGN/star 
formation durations). Old AGNs (e.g., in the last quarter of the durations) that have spent 
long time since the ignition of coevolution would have already consumed the gas through 
the past star formation, with only a tiny amount of leftover for further bulge growth. 

Here, we notice a subsample of AGNs, narrow-line Seyfert 1 galaxy (NLS1), as possibly 
young AGNs. NLSls are characterized (see Pogge 2000) by (i) narrower Balmer lines of 
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hydrogen (FWHM of H/3 < 2000 km s 1 ), relative to usual broad-line Seyfert 1 galaxies 
(BLSls) and QSOs having FWHM of H/3 > 5000 km s" 1 (Osterbrock & Pogge 1985). (ii) 
They often emit strong optical Fe II multiplets (e.g., Halpern & Oke 1987), and show (iii) 
steep and luminous soft X-ray excess (Pounds et al. 1995; Otani et al. 1996; Boiler et 
al. 1996; Wang et al. 1996; Laor et al. 1997; Leighly 1999b). (iv) Rapid X-ray variability 
is another characteristic of NLSls (Otani et al. 1996; Boiler et al. 1997; Leighly 1999a; 
Hayashida 2000). Characteristics in (ii)-(iv) indicate that NLSls are non-obscured objects, 
alike type 2 (i.e., obscured) AGNs. Currently, the most promising picture of NLSls is that 
they contain relatively less massive black holes (with M B h ~ 1O 6_8 M ) and higher L/L E dd 
ratio (~ 1; e.g., Brandt & Boiler 1998; Hayashida 2000; Mineshige et al. 2000), where L 
and LEdd are the AGN luminosity and Eddington luminosity, respectively. NLSls occupy 
10-30% among type 1 AGNs in the local universe (see Kawaguchi et al. 2004b). 

High L/Z/Edd ratios of NLSls imply that their BHs are now growing rapidly via super- 
Eddington accretion (Rees 1992; Mathur 2000; Kawaguchi et al. 2004b). Provided that a 
gas accretion rate onto the central BH, M, during an AGN phase does not drastically change 
with time (see Collin & Kawaguchi 2004), large M/Mbh ratios of NLSls also mean that the 
elapsed time since the ignition of an AGN activity is shorter than that of BLSls. (If the gas 
accretion lasts long enough, then M B h would not be so small.) Bulge luminosity of NLSls 
are smaller than those of BLSls in general (Botte et al. 2004). 

NLSls are likely young active galactic nuclei [during the first 10-30% of the whole 
AGN lifetime (~ 10 8 yr; see §6)], in the process of coevolution of supermassive BHs and 
bulges. This is our working hypothesis in this study. We consider that they are the best and 
unique targets to explore how the BH-bulge coevolution proceeds. If the growth phases of 
BHs and bulges are really overlapped, bulges of NLSls may also be rapidly growing with a 
significant amount of molecular gas accumulated in their bulges. Do they show large ratios 
of molecular-gas to bulge masses, compared with BLSls (which are expected to be older 
than NLSls, based on their lower M/Mbh ratios)? To examine these conjecture, we need 
copious data for numerous NLSls, collected under the condition that the spatial resolution 
is less than the bulge size (typically in the order of lkpc). 

However, NLSls are rare among Seyfert 1 galaxies, and thus we do not have many 
NLSls, whose CO emission are detected, in the nearby universe. Therefore, observations to 
increase the number of CO-detected NLSls are fundamental for listing strong CO emitters, 
with an eye to follow-up sub-kpc observations on them. The MeH^bulge mass relation 
of NLSls and its comparison with the relation of inactive galaxies are subjects of debate 
(Komossa & Xu 2007, and references therein). Our observations would also help to clarify 
if NLSls are really young AGNs. 
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In this paper, we present the observational results of 12 CO (1 — > 0) line from 2 NLSls 
with 5-6" (2-4 kpc at their distance) spatial resolution using the Nobeyama Millimeter Array 
(NMA) at the Nobeyama Radio Observatory. The next section outlines the targets selection, 
and §3 describes the log and data reduction of our radio observations. We then show the 
results in §4, followed by discussions on bulge-mass to stellar-mass ratios in §5. The final 
section is devoted to discussions and summary. Cosmological parameters adopted in this 
paper are H Q = 70kms _1 Mpc -1 , Qm = 0.3 and fl^ = 0.7. 



2. Selection of Targets 

So far, there are only a few NLSls whose CO lines are detected in the nearby uni- 
verse. Intending to increase the number of such NLSls, we thus carried out the following 
observations (with beam size ~6") onto two NLSls with no CO observations reported. 

In order to select targets, we investigated flux densities at 100/im (measured with IRAS 
satellite), which roughly trace the CO line flux (Solomon et al. 1997), for all NLSls at Decl. 
> —25° and redshift z < 0.067 (the sample summarized by Ohta et al. 2007). 

For NGC 4051, of the largest 100/im flux (24 Jy), its CO emission line is clearly detected 
(Young et al. 1995; Maiolino et al. 1997; Kohno et al. 1999a, b). Whereas only an upper- 
limit is obtained for the second 100/im brightest (4.7 Jy) object (Mrk 766 = NGC 4253), 
with M(H 2 ) < 10 7 ' 7 M Q (Taniguchi et al. 1990). It then turned out that the third and forth 
brightest NLSls (~4 Jy) have no radio observations made on their molecular gas. We select 
these two SBb galaxies, IRAS 04312+4008 (z =0.020) and IRAS 05262+4432 (z =0.032), to 
observe. 

We regard these two NLSls as the best targets to elucidate how bulge growth is accom- 
panied with BH growth: By assuming a viliarized Broad Line Region (Wandel et al. 1999; 
Kaspi et al. 2000), M B h is estimated with H/3 width: 



M RH = 4.82 x 10 6 M m 



FWHM(H/?) 
1000 km s" 



1 2 



vL v (B100k) 
10 44 erg s- 1 



0.7 



Although some improved equations in place of the equation above are suggested by, e.g., 
Onken et al. (2004), Kaspi et al. (2005) and Collin et al. (2006), we adopt eq. [I] in this paper 
for simplicity. The FWHM(H/3), optical luminosity and estimated BH masses for the two 
objects are (860 kms" 1 , 10 44 - 6 ergs" 1 , 1O 7 -°M ) and (740 kms" 1 , 10 45 - 2 ergs" 1 , 1O 7 - 3 M ) for 
IRAS 04312 and IRAS 05262, respectively (Ohta et al. 2007). Based on the M BH -estimations 
in this way and our spectral models for M-evaluation (Kawaguchi 2003), they have extremely 
small M-q-r/M ratios (< 2 Myr) even among NLSls [all CO data for NLSls taken from our 
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observations and from literatures will be given in the forthcoming paper (Kawaguchi et al. 
in prep.)]. Thus, they probably have spent little time since each BH began to grow. 



In April-March 2006, aperture synthesis observations toward 2 NLSls were carried out 
with the NMA D-configuration, with the longest baseline of 82m. The beam size (~ 5" x 6") 
corresponds to 2-4 kpc at their distances (see §5 for their bulge effective radii). The NMA 
consists of six 10m antennas, although one antenna was sometimes out of order during our 
observations. 

In all observations, 12 CO (1 — ► 0) (A rcst = 2.6008 mm or z/ rest = 115.271 GHz) lines were 
searched for, by setting the central frequency of each observation to the redshifted frequency. 
The backend we used was the Ultra- Wide-Band Correlator (UWBC; Okumura et al. 2000), 
which was set to cover 512MHz (~1300 kms" 1 ) with 256 channels at 2 MHz resolution (with 
~10km/s effective resolution after the Hanning window function applied). The quasar 
3C111 was used for phase and amplitude reference calibrations, while 3C454.3 was used for 
bandpass calibration. Total on-source integration time on each object and beam size are 
summarized in Table 1. To calibrate absolute flux-scale, Uranus was used. The uncertainty 
in the absolute flux scale is typically ~ 10%. Standard data reduction was performed using 
the UVPROC-II software package developed at NRO (Tsutsumi et al. 1997). Some parts 
of data taken under a bad radio seeing condition were removed from our analysis. After 
clipping a small fraction of data with unusually high-amplitude, the data were then imaged 
with natural UV weighting, and CLEANed, with the NRAO AIPS package. The absolute 
positional accuracy is less than ~ 5 % of the beam size. 



The CO emission from the 2 IRAS galaxies was detected, for the first time, by our 
observations. 



3. 



Observation and Data Reduction 



4. 



Results 



4.1. 



Distribution of Gas 



To show the distribution of the molecular gas of the two targets, we here show the 
CLEANed maps binned over broad velocity ranges (320 or 190 km s _1 ) in figuresHHSJ Through- 
out, north is at the top and east to the left in maps. The ranges integrated are shown in 
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figures [3HU with labels "map". In both objects, the peak flux densities amount up to ~ 16a 
(~ 100 mJy per beam). By looking at the 2MASS images, we find that the peak of CO emis- 
sion from each galaxy coincides with the peak position of the NIR emission of each galaxy. 
The peak position of the CO emission is marked with the crosses, to refer the central position 
in each channel map later (figures [5H6]) . Optical images taken at University-of-Hawaii 88 
inch telescope at J-band, taken under ~ 0.6" (FWHM) image quality (Ohta et al. 2007), are 
also shown to present the morphology of the host galaxies. 

The gas distribution of IRAS 04312 is quite centrally concentrated (fig. [p. The FWHM 
of this signal is about 7.2" x 6.8" . Considering the beam size (5.8" x 5.2"), the spatial extent 
(in diameter) of the gas is estimated to be ~ 4.3" x 4.4" (~ 1.7-1.8 kpc), well smaller than 
the beam. On the other hand, the molecular gas of IRAS 05262 (fig. |2j) seems extended larger 
than the beam size (6.1" x 5.2"). The FWHM of its central CO emission, ~8.3"xl0.7", 
suggests the actual size to be ~5.6"x9.4" (3.5x5.9 kpc), elongated in the North-South 
direction. In §5, these sizes are compared with their bulge effective radii, estimated by a 
couple of methods, to discuss how much fraction of the detected CO gas are located within 
the bulge. 

Incidentally, an isolated CO emission is seen ~ 7" (~ 4.4 kpc in projected distance) south 
from the center of IRAS 05262, which seems to arise from a star forming region evident in 
the optical image. 

4.2. CO Spectra and Molecular-Gas Masses 

Next, we present spectra of CO (1 -> 0) line for the 2 NLSls, in figures [3] (IRAS 04312) 
and H] (IRAS 05262). Spectrum of each central region shows a double-horned shape, indicat- 
ing a rotating molecular disk in origin. The horizontal solid lines, labeled as "map", show 
the ranges used in figures [TH2] to draw the frequency-integrated maps. Downward arrows 
indicate the recession velocity of the targets. 

The CO-to-H2 conversion factor is estimated to be 0.8 and 2.9 M (K km s _1 pc 2 ) -1 
for ultra-luminous infrared galaxies (Downes & Solomon 1998) and for our Galaxy (Dame 
et al. 2001), respectively. When converting CO line flux to H 2 mass, M(H 2 ), the conversion 
factor = 1 M (K km s" 1 pc 2 ) -1 is assumed throughout this paper. The same conversion 
factor is applied for re-computations of M(H 2 ) for objects taken from literatures. 

The velocity-integrated flux from IRAS 04312 measured over the same area used to 
draw figure El is 57±3 Jy kms -1 , meaning that M(H 2 ) = (1.0 ±0.05) x 1O 9 M . The velocity 
integration range is shown in figure [3] with a label "flux". While, the velocity-integrated 
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flux from the central region of IRAS 05262, measured in the box including the CO peak 
(figfJD, is 62±3 Jy kms" 1 , corresponding to M(H 2 ) of (2.8 ±0.1) x 10 9 M Q . These molecular- 
gas masses are the 2nd and 4th largest masses (to our best knowledge) among NLSls, 
with IZwl [M(H 2 ) = 1O 95 M and z = 0.061; Maiolino et al. 1997] and PG 1440+356 
\M(H 2 ) = 1O 9,3 M and z = 0.079; Evans et al. 2001] being the 1st and 3rd ones. Summary of 
all NLSls whose CO emission are detected will be given in the forthcoming paper (Kawaguchi 
et al. in prep.). If we restrict ourselves at z < 0.05 (where 1" angular resolution is equivalent 
to lkpc spatial resolution), these masses are the largest ones among NLSls. 

In the bottom panel of figure HI the spectrum of the isolated CO emission (~ 7" south 
from the center) is also shown. This quite narrow emission line has a velocity-integrated 
flux of 21±1.4 JyKrns" 1 , corresponding to M(H 2 ) = (9.3 ± 0.6) x 1O 8 M . Translating this 
H 2 mass into a star formation rate (SFR) by using equations [5] and [6] shown later, we get an 
SFR of 3.0 Moyr- 1 . 



4.3. Kinematics 

Figures [5] and [6] show channel maps of the two IRAS galaxies. Each map presents a 
map binned over 42kms _1 , and velocity step between adjacent maps is 21 kins" 1 . The cross 
drawn in each panel represents the peak position of the full CO emission (Figs. [Hand [2]). In 
IRAS 05262, a molecular-gas clump at ~ 7" south from the center is evident. 

In both objects, rotating dynamics around each center (cross) is clearly seen. Position- 
velocity diagrams shown in figures [7HS] also indicate rotation. The position angle used for 
each diagram is determined as follows. Two maps (of high S /N peaks) are chosen from each 
channel map. Then, the peak positions in the two maps (with two channels separated by 
~ 210 km s -1 ) are connected by a line. 

Now, we derive dynamical mass (i.e., enclosed mass) Md yn from these position-velocity 
diagrams. We estimate an inclination angle i for each galaxy by measuring the ellipticity 
of the /-band image (Ohta et al. 2007), where the inclination angle of deg means face-on 
view on a rotating CO gas disk: i = 40° ± 1° (IRAS 04312) and i = 53° ± 5° (IRAS 05262). 
Then, the dynamical mass is evaluated as 

M dr ^5.8xl0»M 9 ( I5 ^- T ) 2 ( i |-)( sml ,-, (2) 

where AV and R are the velocity difference and the separation. Positional separations R and 
velocity differences AV between two peaks in the position-velocity diagrams are 2.7" (1.1 kpc 
in projected distance) and 164 km s -1 for IRAS 04312, and 5.1" (3.2 kpc) and 165 km s -1 for 
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IRAS 05262, respectively. The enclosed mass Md yn within the central region of IRAS 04312 
is 4.1 x 1O 9 M . While for IRAS 05262, M dyn = 8.0 x 1O 9 M . The peak-to-peak separations 
represent the lower limits of their CO-gas sizes. As we saw in § 4.1, the CO gas in IRAS 04312 
is indeed more concentrated than that in IRAS 05262. 



We here discuss the sizes of the CO gas in two NLSls (§4.1 and §4.3), and compare 
them with their bulge sizes. The spatial extent (in diameter) of the CO gas in IRAS 04312 is 
estimated to be ~4.3" x4.4" (~ 1.7-1.8 kpc). On the other hand, IRAS 05262 seems to have 
an extended and elongated CO gas at its center with the size of ~5.6" x9.4" (3.5x5.9kpc). 
Peak-to-peak separations in position-velocity diagrams (Figs [7] and [8]), 2.7" (=1.1 kpc) and 
5.1" (=3.2 kpc) for IRAS 04312 and IRAS 05262, respectively, also indicate that the molec- 
ular gas of IRAS 05262 is more extended than that of IRAS 04312. 

Now, let us estimate the bulge effective radius (Reg) from [O III] width as follows. The 
FWHM of [O III] emission lines of IRAS 04312 is 380 km s" 1 (Veron-Cetty et al. 2001). Since 
[O III] widths and velocity dispersions (a*) are nearly the same (Heckman et al. 1989; Nelson 
& Whittle 1996), a* for this object is evaluated to be about 160 km s -1 , by dividing the [O 
III] FWHM width by 2.35. Elliptical galaxies and bulges of disk galaxies obey a canonical 
relation about their dynamics, size, etc (the fundamental plane). The velocity dispersion 
can be translated into the size (Guzman et al. 1993): 



Thus, the bulge effective radius of IRAS 04312 is estimated to be about 2.6 kpc or 6.4" at its 
distance. However, the inferred bulge size (13" in diameter) is comparable to the length of 
the bar seen in the optical image in the North-South direction. Therefore, this method likely 
overestimates R e g for this object. Adopting the same method to IRAS 05262 by taking its 
[O III] width (365 km s" 1 ) from Veron-Cetty et al. (2001), its R e g is about 2.3 kpc or 3.6" 
at its distance. 

We note that there is a large scatter about the regression equation (eq. [3]). Thus, R e s 
estimation using eq. [3] alone should not be taken too much seriously. For instance, our visual 
investigations on the radial profile of the two IRAS galaxies (data presented by Ohta et al. 
2007) lead us rough estimations of their R e g to be ~ 0.7" and 3.1" for IRAS 04312 and 
IRAS 05262, respectively. Because of strong emission from the nuclei, a large uncertainty of 
R e g evaluated by this method is also inevitable. 



5. Sizes of CO Gas and Bulge 




(3) 
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For IRAS 04312, two methods above provide so different estimations (13" and 1.4" in 
diameter). Thus, it is premature to compare these with the size of the CO gas (~ 4" x 4"), 
and to discuss how much fraction of the total CO flux comes from its bulge region. While, 
two values are consistent each other for IRAS 05262 (7" and 6" in diameter). Since the 
spatial extent of the detected molecular gas in IRAS 05262 (~5.6" x9.4") seems a bit larger 
than the estimated bulge size, the molecular gas harbored within the bugle would be slightly 
smaller than the detected M(H 2 ). To the first order, the observed M(H 2 ) well represents 
the M(H 2 ) of the bulge of IRAS 05262. It is unclear whether stars newly born from the 
molecular gas belong to the bulge or to the galactic disk. In this study, we regard star 
formation within R e g as bulge growth. 



6. Molecular Gas and Bulge Masses 

Here, we compare the molecular-gas masses of these two IRAS galaxies with the dy- 
namical masses (M dyn ) and with the bulge masses (M bulge ). These masses are summarized 
in Table 2. 

The dynamical mass Md yn is inferred from the CO-gas kinematics. We discuss the 
position-velocity diagrams (figures [7] and ED • The enclosed mass M<j yn within the central 
region of IRAS 04312 is 4.1 x 1O 9 M , resulting in a M(H 2 )/M dyn ratio of 0.25. While for 
IRAS 05262, M dyn = 8.0 x 1O 9 M , and the M(H 2 )/M dyn ratio is then 0.35. If we estimate 
the current stellar mass as Md yn — M(H 2 ), the molecular gas masses for these two objects 
account to 33 % and 54 % of the current stellar massed- However, all the molecular gas 
will not be consumed by star formation during this star formation episode (in ~ 10 7 ~ 8 yr), 
since a molecular cloud is to be destroyed by photoevapolation, mass loss and supernova 
explosion (Williams & McKee 1997). Recycling of the remaining diffuse gas may occur, but 
on a timescale much longer than the timescale we are considering. Since we are interested 
in the BH-bulge evolution during one coevolution episode, we dot not go further upon the 
recycling process (on a long timescale) in this study. The fraction of a molecular gas mass 
that is converted into stars in molecular clouds (efficiency of star formation) ranges from a 
few percent (Myers et al. 1986; Wilson & Matthews 1995) to tens percent (Hillenbrand & 
Hartmann 1998; Lada & Lada 2003). Adopting a star formation efficiency of 0.1, new stars 
with masses of 3 — 5% of the preexisting stars will be formed in the near future (e.g., within 
the timescale of cloud-destruction, ~ 10 7 5 yr; Williams & McKee 1997) from the molecular 



1 HI gas is normally negligible (compared to H2 masses) within several kpc ("molecular front") around 
galactic centers (Sofue 1994; Sofue & Nakai 1994). 
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gas. 



As we discussed in §5, the bulge size is not definitively determined. Thus, we can not be 
completely sure that the dynamical mass (or virial mass) equals to the bulge mass. Therefore, 
we next try to have another estimation of the bulge mass as a complementary measure. 
The bulge mass, Mb u i ge , can be estimated from kinematics of stars/clouds embedded in 
the gravitational potential of the bulge, in the form where M bulge oc (line width) 4 . To 
determine the normalization in this relationship, we refer to the M — a* relationship found 
in normal galaxies. If we adopt the two relations, Mbh = 8 x 10 6 M q [^/(lOOkms -1 )] 4 and 
M bu ig e = 710 M BH (e.g., Tremaine et al. 2002; Gebhardt et al. 2000; Haring & Rix 2004), we 
obtain the following relationship: 



The scatter in Mbh^o"* relation is about 0.3 dex (Gebhardt et al. 2000), whereas the uncer- 
tainty in the Mb u i ge /MBH ratio is also ~ 0.3 dex (e.g., Marconi & Hunt 2003). The uncertainty 
of the bulge masses evaluated in this way would thus be about 0.4 dex. 

CO widths might be useful to evaluate ex* and then bulge masses (Shields et al. 2006; 
Ho 2007; see § IT. 2j) . We here follow their procedure. Assuming that the dispersion of 
CO emission line exco equals to ex*, as suggested by Shields et al., we adopt aco of 109 
and HOkms" 1 for IRAS 04312 and IRAS 05262, respectively (see gZ2J). With eq. HJ bulge 
masses are estimated as M bu i ge = 8.0 x 1O 9 M and 8.4 x 1O 9 M , resulting in M(H 2 )/M hu \ gc 
ratios of 0.13 and 0.33, respectively. [If [O III] widths are used as surrogate for cx* in eq. HJ 
M{H%)j Mb u ige ratios decrease by a factor of 4 (see § 7.2).] Star formation will undergo with 
1 — 3% of the masses of the current stellar masses for the star formation efficiency of 0.1. 
Even if we again estimate the mass of preexisting stars by M bu ig e — M(H 2 ) to account for 
the contribution of gas mass, these numbers do not change dramatically. 

Let us now summarize this section. Mass of new stars that will be formed (within 
~ 10 7 ' 5 yr from now on) from the molecular gas is not a huge fraction of the current stellar 
masses (1-5%). Here, we conclude that the bulges of these two NLSls do not evolve so 
much in this timescale. If newly born stars are to belong to the galactic disks rather than 
the bulges, the conclusion (small growth in the bulge mass) is enhanced. 

This result can be interpreted in some ways. (1) Coevolution is not simultaneous at all: 
bulges (and elliptical galaxies) were formed in the early universe (Akiyama 2005; Woo et al. 
2006), while BH growth via accretion still happens today. (2) Alternatively, coevolution 
of BHs and bulges can be (quasi-) simultaneous with two growth episodes overlapped or 
close in time. (2-1) Our result indicates that the bulge growth have finished prior to the 
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BH growth (e.g., Umemura 2001). (2-2) Or, coevolution may occur intermittently through 
a number of episodic events. Each coevolving AGN could have very short durations of BH 
and bulge growth, with small increases in Mb u i ge and in Mbh, resulting in a small fraction of 
molecular gas associated with the AGN even in the early evolutionary stage (NLSls). The 
total lifetime of an AGN is estimated as ~ 10 8 yr (Martini 2004; Jakobsen et al. 2003). 
While, a lower limit for the duration of an episodic event of an AGN is about 10 4 yr, based 
on the proximity effect (e.g., Bajtlik et al. 1988) and the sizes of ionization-bounded narrow- 
line regions (Bennert et al. 2002). Therefore, BH-bulge coevolution could be divided into 
numerous shorter events. 



■>FIR — 19 Lq 



(5) 



7. Discussions and Summary 

7.1. Star Formation Rates v.s. Accretion Rates 

Here, we discuss the relation between star formation rates and accretion rates, and the 
^buige — -Mbh relation. 

From FIR luminosity or from molecular-gas mass, the current star formation rate can 
be estimated. The regression relation between 40-120/im FIR luminosit}@ and M(H 2 ) is 
typically (for galaxies similar to the Milky Way) expressed as (see Kennicutt 1998; Tinney 
et al. 1990), 

~ M(H 2 ) ' 

M Q 

The scatter around this linear relation is about ±0.3 dex (Tinney et al. 1990). (Again, the 
regression is re- normalized here by adopting the same conversion factor as we used so far.) 
Now, we estimate the FIR luminosity at each central region (free from the contribution of 
AGNs) for the two NLSls, based on their molecular gas masses. (On the other hand, observed 
FIR luminosity by IRAS, shown below, includes the FIR emission from entire galaxy and 
from AGN.) Inferred masses (§4.2) indicate that their Lfir are 1.9 x 10 10 L & and 5.3 x 10 10 
L Q , respectively. 

By translating these FIR luminosities to the star formation rate (SFR) by (Kennicutt 
1998) 

(6) 



SFR= l.TMoyr- 1 ( -^-^ 
oy \10 10 L e ) 



2 Lfir = 0.65 x uL v (60^m) + 0.42 x uL u (100[im) (Hclou ct al. 1985). 
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their SFRs (at their central region) are expected to be 3.3 and 9.2 M Q yr _1 , respectively. 
Most of normalizations for eq. [6] in literatures lie within ± 30 % for starbursts. Disk galaxies, 
on average, may have a normalization a factor of 1.8 higher than eq. [6] (Buat & Xu 1996). 
Taking all the uncertainties of normalizations in eqs. [5] and El an error of SFR estimated 
from M(H 2 ) would be ±0.3 + 0.3 dex. 

Incidentally, observed Lfir of the 2 IRAS galaxies (measured via IRAS, including con- 
tributions from entire host galaxy and from nuclei) are 2.7 x 10 10 L Q (IRAS 04312) and 6.3 
x 10 10 L Q (IRAS 05262). If these L F ir are used instead, their SFRs are evaluated to be 
5 and 11 M yr _1 , respectively. Since these values likely overestimate the central Lfir and 
thus central SFRs, we hereafter refer the SFRs estimated in the previous paragraph. 

Next, we evaluate the accretion rates onto their central BHs M , by adopting two simpli- 
fied methods. Firstly, we estimate M based on bolometric luminosities, with an assumption 
of a constant radiation efficiency Bolometric luminosity is assumed to be 10 times the op- 
tical luminosity, i>L u (hlWA) (e.g., Elvis et al. 1994). [We note that up to ~ 50% of the 
optical luminosity may originate in the host galaxies (Surace et al. 2001; Bentz et al. 2006). 
Thus, discussions below overestimate M by up to a factor of ~ 2.] The radiation efficiency 
(with respect to the rest-mass energy of the infalling gas) is set to be 0.05 (e.g., Shakura & 
Sunyaev 1973; Novikov & Thorne 1973). Then, M are estimated to be 1.3 and 5.2M /yr, 
resulting in the SFR/M ratios of 2.5 and 1.8 for IRAS 04312 and IRAS 05262, respectively. 

When M exceeds the Eddington limit (~ 16 L^ddfc 2 ), however, the radiation efficiency 
drops due to the onset of photon trapping (Begelman 1978; Abramowicz et al. 1988). Thus, 
estimates using accretion disk models are more adequate than those using a constant effi- 
ciency (see e.g., Kawaguchi 2003; Collin & Kawaguchi 2004). Here, M is then estimated 
based on optical luminosity and the (standard) accretion disk model with a face-on disk 
geometry assumed (e.g., Bechtold et al. 1987)0: 

M 



M yr 



-i 



2.1 



( M] 



BH 



Vl0 7 M Q 



^(51001) 
10 44 erg s- 1 



1.5 



(7) 



3 The photon trapping effect takes place only in the inner part of the disk where UV and soft X-ray 
photons emerge. In outer region of the disk which radiates optical continuum emission, self-gravity (rather 
than the gravity from the central BH) governs the disk dynamics (Kawaguchi et al. 2004a). Although the 
standard accretion disk model (Shakura & Sunyaev 1973) works only inside the self-gravity radius, we here 
adopt the standard disk model beyond the self-gravity radius for the sake of its simplicity. 
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With this equation, we get 19 and 76 MqJt' 1 for IRAS 04312 and IRAS 05262, respectively. 
Here, we obtain yet smaller SFR/M ratios of 0.17 and 0.12, respectively. 

Then, let us discuss these ratios in the context of Mt, u i ge — Mbh relation. The ratio of 
bulge masses to BH masses in galaxies is about 700 (Haring & Rix 2004). This number is 
considerably larger than the SFR-to-M ratios for the two NLSls (~ 1) estimated above. 
What does this huge difference mean then? 

Our result (the SFR-to-M ratios ~ 1 for the two NLSls) can indicate that M in NLSls 
is temporarily large. To compensate the intense BH growth during the NLSls, a short-term, 
intense starburst (with much higher SFR than the current SFR) happens before or in the 
future (i.e., bulge growth and BH growth phases appear without an overlapped period). If 
the duration of star formation is much longer than the duration of AGN phase instead, such 
an intense starburst is not required. 



7.2. CO Width as Surrogate for aj 

Heckman et al. (1989) showed that CO and [O III] lines have similar widths. While, 
[O III] line widths are now often used as surrogate for bulge stellar dispersion a* of AGNs 
(Nelson & Whittle 1996; Nelson 2000; Shields et al. 2003; Boroson 2003). Then, Shields 
et al. (2006) proposed CO width as surrogate for ex* (see also Ho 2007), although they found 
a systematic difference (cr[ iii]/cr C o ~ 0.15 dex, on average, for the seven PG quasars they 
looked at). By using the CO width for high redshift quasars, they examined the evolution of 
M — a* relation at z > 3, arguing that high-z quasars have smaller M bulge /M B H ratios than 
quasars at z < 3. 

Here, we compare CO width and [O III] width using our results, in order to assess 
the relation between CO width and [O III] width for NLSls. Full- Width at Zero-Intensity 
(FWZI), evaluated with the ranges drawn in figures [SHU with "flux" labels, for the two 
NLSls are 384 km s" 1 (IRAS 04312) and 388 km s" 1 (IRAS 05262). Following Shields et al. 
(2006), we also assume that FWHM(CO) « § FWZI. Then, FWHM(CO) for these two IRAS 
galaxies are about 260 km s -1 , which is less than their FWHMQOIII]) widths (380 km s" 1 and 
365 km s -1 , respectively). This is the same trend that Shields et al. (2006) presented for PG 
quasars. Therefore, usage of CO width to estimate a* can lead a systematic error upon Mb u i ge 
for NLSls as well. The difference between the CO and [O III] widths are 0.15 — 0.17 dex for 
the two NLSls, providing a factor of 4 less M bu i ge for CO-based mass. 
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7.3. Summary 

The tight correlation between mass of central BHs and mass of galactic spheroid compo- 
nents (galactic bulge or elliptical galaxy itself) in galaxies and AGNs indicates contempora- 
neous coevolution of the two massive systems. However, no direct evidence for simultaneous 
growth has been found. We aim to reveal whether BH growth and bulge growth take place 
(quasi-) simultaneously, by determining the amount and distribution of molecular gas (fuel 
for the coming star formation) in young AGNs. 

NLSls are characterized by relatively small BH masses and high Eddington ratios. 
Larger M/M B h ratios of NLSls than those of BLSls indicate that BHs in NLSls are now 
growing rapidly, and that NLSls are younger than BLSls. If the growth phases of BHs and 
bulges are really overlapped, bulges of NLSls may also be rapidly growing with a significant 
amount of molecular gas accumulated in their bulges. However, NLSls are only 10-30% 
among Seyfert 1 galaxies, and thus we do not have many NLSls, whose CO emission is 
detected, in the nearby universe. In order to investigate the nature of molecular gas in 
NLSls in general, it is fundamental to increase the number of nearby NLSls with their CO 
emission lines detected. 

Thus, we performed 12 CO (1 — > 0) observations using the NMA D-configuration for 
two FIR-bright NLSls located at z = 0.020 - 0.032 with 5-6" (2-4 kpc at their distance) 
resolution. We detect their CO emission for the first time, with M(H2) of (1 — 3) x 1O 9 M . 
These molecular-gas masses are the 2nd and 4th largest masses among NLSls. 

By estimating M dyn and Mb u i gc with the CO channel maps and CO line widths, we 
get M(H 2 )/M dyn and M(H 2 )/M hu i ge ratios of 0.25-0.35 and 0.13-0.33, respectively. Taking 
account the star formation efficiency (~ 0.1), the increase in Mb u i ge in those NLSls in the 
near future is expected not to be a huge fraction (1-5% of the mass of the preexisting 
stars). This result can be interpreted in some ways. Coevolution may not be simultaneous 
at all: bulges (and elliptical galaxies) were formed in the early universe, while BH growth 
via accretion still happens today. Or, the bulge growth may have finished prior to the BH 
growth. Alternatively, coevolution may proceed with many, occasional discrete events, where 
one coevolution event produces only a small amount of mass growth of BHs and of bulges. 

We further discuss the relation between star formation rates and accretion rates, in the 
context of the MBH _ -^buige correlation. The two NLSls turned out to have much smaller 
SFR/M ratios (~1) than the Mb u i ge /^BH ratios of normal and active galaxies. This huge 
difference can indicate either that bulge growth and BH growth phases appear without an 
overlapped period, or that the duration of star formation is much larger than the duration 
of AGN phase. 
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For the two NLSls, CO widths are less than their [O III] widths (by 0.15-0.17 dex), 
which is the same trend that Shields et al. (2006) presented for several PG quasars. Therefore, 
usage of CO width to estimate o* can lead a systematic error (a factor of 4 less masses) upon 
Mauigc for NLSls as well. 

We thank the NMA staff for helping observations. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory California Institute of Technology under contract with the National Aeronautics 
and Space Administration. T.K. thanks the financial supports from a Grant-in-Aid for 
Scientific Research of JSPS (18840038) and from The Research Institute of Aoyama Gakuin 
University. The Nobeyama Radio Observatory is a branch of the National Astronomical 
Observatory of Japan, the National Institutes of Natural Sciences (NINS). 
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Table 1. Log of Observations 



Target 


rcdshift 


Integration Time 


Beam Size 


Scale 






[hours] 


["x"] 


[kpc/1"] 


IRAS 04312 . . . 


0.020 


6.6 


5.8 x 5.2 


0.40 


IRAS 05262 . . . 


0.032 


7.7 


6.1 x 5.2 


0.63 
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Table 2. Summary of Various Masses 



Object 


M{H 2 ) 


-Mdyn 


Ratio a 


Mbulgo 


Ratio" 




[10 9 M Q ] 


[1O 9 M ] 




[1O 9 M ] 




IRAS 04312 . . . 


1.0 


4.1 


0.25 


8.0 


0.13 


IRAS 05262 . . . 


2.8 


8.0 


0.35 


8.4 


0.33 



a Ratio of the molecular gas masses MiH^) with respect to the dynam- 
ical masses Md yn or bulge masses Mb u i g e- 
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Fig. 1. — CLEANed contour map of CO (1 — > 0) emission line from IRAS 04312 binned in 
frequency over 320 km s -1 , on top of the optical image (Ohta et al. 2007). Contour levels 
show —2, 2, 4, . . . 14, 16 a, with la = 2.0 Jy beam -1 km s -1 . The 4" x4" cross represents the 
peak position of the CO distribution. At the left-bottom corner, the beam size (5.8" x 5.2") 
is shown. 
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Fig. 2. — The same as figure [TJ, but for IRAS 05262. Larger square including the galactic 
center and another square, south from the center, indicate the areas over which the spectra 
(figureHJ) are integrated. The contour map shows the CO emission integrated over 190 km s -1 , 
with la = 1.2 Jybeam -1 kms -1 . The beam size is 6.1" x 5.2". 
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Fig. 3. — The CO (1 — > 0) spectrum of the central region of IRAS 04312, integrated over 
12"xl2" box centered on the peak of the CO emission. (Hereafter, velocity notation is 
based on optical convention.) The bin-width and frequency-step are the same as those for 
the channel maps (figJSJ). Downward arrow indicates the recession velocity of the galaxy, 
measured with H I gas (Paturel et al. 2003). The horizontal solid line indicates the range 
used in figured] to draw the frequency-integrated map. Another horizontal line with a label 
"flux" indicates the range for measuring the velocity-integrated flux. 
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Fig. 4. — The same as figure [3j but for IRAS 05262. Top panel shows the spectrum of the 
central CO gas, while the bottom presents the spectrum of the H II region ~ 7" south from 
the center. 
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Fig. 5. — CLEANed channel map of CO (1 — > 0) emission line around IRAS 04312 over 
30" x 30". Each map is binned in frequency over 42kms _1 , showing (-3, -2, 2, 3, 4, 5 
. . .) x cr with lor = 19 mJy per beam, and drawn every 21kms _1 step. Crosses represent 
the peak position of the full CO emission (Fig. [1]). Number at upper-right corner of each 
panel represents the velocity offset (in kms" 1 ) from the recession velocity of the galaxy 
(6145 km s" 1 ). 




Fig. 6. — The same as figure[5l but for IRAS 05262 (with a recession velocity of 
Here, la = 14mJy per beam. 



9646 km s" 1 ). 



-28- 



6350 



6300 



6250 



6200 



6150 



u 
_o 

cu 

> 

1/1 



6100 



6050 



6000 



5950 




East 



West 



-10 



-5 5 

Angular Distance (arcsec) 



10 



Fig. 7. — Posit ion- Velocity diagram of 12 CO (1 — > 0) emission from IRAS 04312 along the 
position angle of 108°. The abscissa is the angular distance from the CO emission peak. The 
contour levels drawn are the same as those of the channel map for the same object (fig 15]). 
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Fig. 8. — The same as figured but for IRAS 05262 (with the position angle of 97°). 



